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REMARKS 

This document is submitted in response to the Office Action dated December 18, 2006 
("Office Action"). 

Applicants have amended claim 20, to replace "adverse reactions" with "Steven-Johnson 
syndrome or toxic epidermal necrolysis in response to carbamazepine." Support for this 
amendment can be found at, e.g., in original claims 2 and 5. Applicants have also amended 
claims 1 and 20 to remove non-elected subject matter. Further, Applicants have added new 
claim 26, support for which appears, for example, in Table 1, page 30.^ Finally, Applicants have 
cancelled claims 3, 6, and 21. 

Upon entry of the amendments, claims 1, 8-12, 20, and 22-26 will be under examination. 
Applicants respectfully request that the Examiner reconsider this application in view of the 
following remarks. 

Claim Objections 

Claims 1, 12, and 20 are objected to for reciting subject matter that is not elected. See the 
Office Action, pages 2-3, section 1 . Applicants have deleted the non-elected subject matter fi-om 
claims 1 and 20. Applicants, however, disagree that non-elected species should be removed 
from claim 12. Please note that claim 12 recites a Markush group that includes six species of 
equivalent genetic marker of HLA-B*1502, e.g., the elected species HLA-Cw*0801. MPEP § 
803.02 requires that the patentability of the elected species be examined first and if it is found 
allowable, the examination should be extended to the non-elected species. Applicants thus 
submit that it is not necessary to remove the non-elected species fi:om claim 12 at this point. 

Claim Rejection under 35 U.S.C. § 1 12. First Paragraph (Enablement) 

Claims 1, 3, 6, 8-12, and 20-25 stand rejected for lack of enablement on four grounds. 



' Table 1 shows the frequency of HLA-B*1502 in 42 Taiwanese patients with carbamazepine induced Stevens- 
Johnson syndrome/toxic epidermal necrolysis. As Taiwanese are Mongoloids or descendents of Mongoloids, Table 
1 supports new claim 26. 
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Among the rejected claims, claims 3, 6, and 21 have been cancelled. Applicants will traverse 
each of the grounds below. 

I 

Apparently referring to claims 20-25, the Examiner asserts that the Specification does not 
"reasonably provide enablement for assessing the risk of any other adverse reactions in response 
to any other drugs." See the Office Action, page 3, third paragraph; emphases added. Note that 
claim 21 has been cancelled. 

Apphcants have amended claim 20 to limit it to only one adverse drug reaction, i.e., 
Stevens- Johnson Syndrome (SJS)/toxic epidermal necrolysis (TEN), and only one drug, i.e., 
carbamazepone (CBZ). It is submitted that the amendments have rendered moot the Examiner's 
ground for rejecting claim 20, as well as claims 22-25 dependent from it. 

II 

The Examiner asserts that the Specification is not enabling for assessing the risk of 
adverse drug reactions in any human population other than Taiwanese. See the Office Action, 
page 3, third paragraph. Applicants respectfully disagree. 

Independent claim 1 is discussed first. This claim is drawn to a method of assessing the 
risk of a patient for developing SJS/TEN in response to CBZ by detecting the presence of 
HLA-B*1502. 

Citing Lonjou et al, Alfirevic et al., and Hung et al., the Examiner rejects claim 1 on the 
ground that: [i]t is unpredictable as to whether or not HLA-B*1502 would be indicative of risk of 
an adverse drug reaction in another population." See the Office Action, page 6, last paragraph, 
and page 11, second paragraph, through page 12, first paragraph. Applicants disagree. 

Applicants submit that the Examiner has merged two distinct concepts: HLA-B*1502 
being a risk indicator of a patient (required by claim 1) and HLA-B*1502 being a prediction 
marker in a population (discussed in the cited references). As a risk indicator, HLA-B*1 502 is 
used to determine whether a specific patient is at risk for CBZ-induced SJS/TEN. As such, 
HLA-B* 1 502 only needs to be associated with this adverse drug reaction. On the other hand, 
being a prediction marker, HLA-B* 1 502 is used to identify individuals who bear the risk of 
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developing CBZ-induced SJS/TEN in an entire population. Accordingly, the allele needs to be 
present at a relatively high frequency in that population so that it can be detected. As set forth 
below, the cited references merely suggest that HLA-B*1502 may not be a good prediction 
marker to identify risk-bearing individuals in Caucasian population due to its low frequency in 
that population. In fact, none of the cited references discredits using HLA-B*1502 as a risk 
indicator for patients developing CBZ-induced SJS/TEN, as required by claim 1 . 

Lonjou reports that among 12 CBZ-induced SJS/TEN patients (4 Asians and 8 
Caucasians), none of the 8 Caucasian patients has this allele. See Abstract. Accordingly, it 
concludes that "HLA-B*1502 is not a useftil prediction marker of CBZ related SJS/TEN in the 
European population." See page 3, left column, second paragraph; emphasis added. On the 
other hand, it is completely silent on whether a Caucasian patient carrying HLA-B*1502 would 
be at risk for CBZ-induced SJS/TEN, as required by claim 1 . 

Alfirevic reports that (1) none of 56 CBZ-sensitive Caucasian patients (including 2 
SJS/TEN patients) was positive for HLA-B*1502, and (2) none of 43 CBZ-tolerant Caucasian 
patients carry this HLA allele. See Abstract. It thus concludes that HLA-B* 1 502 is not a marker 
for CBZ-induced SJS/TEN in patients of all ethnicities. Similar to Lonjou, Alfirevic only 
addresses whether HLA-B* 1 502 could be used as a prediction marker in an ethnic group in 
which the allele is present at a very low frequency. Again, it does not contradict Applicants' 
position that a patient (in any ethnic group) carrying HLA-B* 1502 is at risk for CBZ-induced 
SJS/TEN. 

Relying on the teachings in Lonjou et al. and Alfirevic, i.e., Caucasian SJS/TEN patients 
being examined do not carry HLA-B* 1502, the Examiner concludes that this HLA allele is not 
associated with the drug reaction in Caucasians. See the Office Action, page 12, second 
paragraph. Applicants would like to point out that HLA-B* 1502 being a risk indicator of CBZ- 
induced SJS/TEN, as required by claim 1, does not compel the conclusion that it is the only risk 
indicator of this drug reaction so that every CBZ-induced SJS/TEN patient carries it. In other 
words, it is likely that Caucasian CBZ-incuded SJS/TEN patients carry risk indicators other than 
HLA-B* 1502. Accordingly, the fact that some CBZ-induced SJS/TEN patients do not carry 
HLA-B* 1502 does not contradict the conclusion that HLA-B* 1502 is a risk indicator of this 
drug reaction, i.e., a patient carrying this HLA allele is at risk for CBZ-induced SJS/TEN. 
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Hung et al. teaches that "[a]n allele associated with a particular phenotype or disease may 
be present in different frequencies in different ethnic groups." See page 233, left column; 
emphasis added. Accordingly, when the frequency of an allele is low in a population, its 
association with a certain disease would be difficult to detect in that population. Indeed, 
according to this reference, "[t]he allele fi-equency of HLA-B*1502 is approximately 4-7% in 
Southern Han Chinese, but only 0-0.1% in Caucasians." See page 233, left column, last 
paragraph through right column, first paragraph; emphases added. This information explains the 
difficulty in detecting the association between HLA-B*1502 and CBZ-induced SJS/TEN in 
Caucasians, which is not contradictory to the teachings of the present application, i.e., a patient 
carrying HLA-B*1502 bears the risk of developing CBZ-induced SJS/TEN. 

As a matter of fact, Hung et al. actually supports the position that HLA-B* 1 502 can be 
used as a risk indicator of CBZ-induced SJS/TEN in any patient. It states that "the HLA-B* 1502 
allele is quite common in Malaysia (approximately 8,4% frequency), and CBZ is reported to be 
the major offending drug for SJS/TEN in that population (35.7% of SJS/TEN patients)." See 
page 233, left colunrn, last paragraph through right column, first paragraph, emphases added. It 
fiirther reports low incidence of CBZ-induced SJS/TEN in Caucasians, who have a low allele 
frequency of HLA-B*1502 (i.e., 0-0,1%). See page 233, left colunrn, last paragraph through 
right column, first paragraph. Based on the above information, Hung et al. proceeds to conclude 
that the HLA-B* 1502 allele frequency positively correlates with the prevalence of CBZ-induced 
SJS/TEN in different populations. See page 233, right colunrn, lines. 8-16. hi other words, 
HLA-B* 1502 highly correlates with CBZ-induced SJS/TEN, and could be used as a risk 
indicator of the drug reaction in any patient. 

In sum, none of the three references relied on by the Examiner suggests that HLA- 
B*1502 is not a risk indicator of CBZ-induced SJS/TEN in a Caucasian patient, i.e., a Caucasian 
patient carrying this HLA allele is not at risk for SJS/TEN. Absent such suggestion, these 
references cannot serve as reasonable bases to challenge the enablement of claim 1 in its entire 
scope, i.e., HLA-B* 1502 can be used to assess the risk of developing CBZ-induced SJS/TEN in 
a patient, regardless of the patient's ethnic backgroxmd. 
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In view of the foregoing remarks, Applicants submit that it is improper to reject claim 1 
for lack of enablement. Nor is it proper to reject claims 8-12, all of which depend from claim 1. 
For the same reasons set forth above, it is also improper to reject claims 20 and 22-25 for lack of 
enablement. 

Ill 

Referring to claim 1 1, the Examiner alleges that "[i]t is unpredictable as to whether or not 
the presence of any equivalent genetic marker is useful for determining the presence of the HLA- 
B*1502 allele or for the assessment of risk of drug adverse reaction." See the Office Action, 
page 8, second paragraph. Applicants disagree. 

The Specification defines an "equivalent genetic marker" of an allele of interest as "a 
genetic marker that is linked to the allele of interest." See page 8, Hues 19-20. It further states 
that "[t]he useful equivalent genetic markers in the present invention display a linkage 
disequilibrium with the allele of interest." See page 8, lines 20-21 . Based on these teachings, 
one skilled person in the art would know that the presence of an equivalent genetic marker of 
HLA-B*1502 indicates the presence of this HLA allele. Thus, he or she would recognize that 
any equivalent genetic marker, as recited in claim 11, could be useful for determining the 
presence of HLA-B* 1 502 and thus for assessing the risk of developing CBZ-induced SJS/TEN. 

For the above reasons, Applicants respectfully request that the Examiner withdraw this 
ground for rejecting claim 1 1 . 

IV 

Finally, referring to claim 12, the Examiner asserts that "there is no statistical analysis of 
the significance of the association of HLA-Cw*0801 with carbamazepine-induced SJS/TEN, nor 
any analysis of the linkage of HLA-B*1502 with HLA-Cw*0801." See the Office Action, 
page 8, second paragraph. Applicants disagree. 

It is well known in the art that HLA-B* 1502 has a strong linkage disequilibrium with 
HLA-Cw*0801. See Romphruk et al., European 1 ImmunogeneL, 30:153-158 (2003), page 4, 
Table 4 and right column, second paragraph; and Romphruk et al., Tissue Antigens, 58:83-89 
(2001), page 86, Table 4, and page 87, left column, first paragraph, (copies submitted herewith 
as Exhibits 1 and 2). Further, the Specification teaches that 38 of 42 carbamazepine-induced 
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SJS/TEN patients (HLA-B*1502 carriers) have HLA-Cw*0801 allele. In other words, as many 
as 90% (38/42) of the SJS/TEN patients carry HLA-Cw*0801 . As the association between 
HLA-Cw*0801 and CBZ-induced SJS/TEN (90%) is clearly significant, one skilled person in 
the art would readily recognize that HLA-Cw*0801, an allele co-present with HLA-B*1502, is 
an indicator of the risk of developing CBZ-induced SJS/TEN in a patient. Indeed, this 
conclusion is supported by a post-filing reference, Hung et al., which shows that 56 out of 60 
SJS/TEN patients carry HLA-Cw*0801 . See Hung et al., Pharmacogenetics and Genomics, 
16(4):297-306 (2006) (copy submitted herewith as Exhibit 3). 

For the reasons set forth above, Applicants submit that the Specification is enabling for 
claim 12, which recites HLA-Cw*0801 . It is thus respectfiiUy requested that the Examiner 
withdraw this lack-of-enablement rejection. 

New Claim 

New claim 26 covers a method of assessing the risk of developing CBZ-induced 
SJS/TEN in a Mongoloid or a descendent of a Mongoloid. The Examiner acknowledges that the 
Specification is enabling for assessing the risk in Taiwanese patients. See the Office Action, 
page 3, third paragraph. As all Taiwanese are Mongoloid, this claim, limited to Mongoloid 
patients, is clearly enabled. 



CONCLUSION 

In view of the above remarks, it is submitted that claims 1,8-12, 20, and 22-26 are 
enabled. Allowance of these claims is respectfiilly solicited. 
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HLA-B""! 5 subtypes in the population of north-eastern Thailand 

^ A. Romphruk/t K. Phongaen,^ J. Chbtechai,^ C. Puapairoj,§ C. Leelayuwat/t & A. V. Romphruktl 



Summary 

The HLA-B*15 group is the most polymorphic HLA-B 
allele and so has several subtypes. These subtypes have 
not been defined in the population of north-eastern Thai- 
land (NET). In a previous study, using polymerase chain 
reaction-sequence-specific primers (PCR-SSP), subtypes 
were categorized into four groups, namely: group I: HLA- 
B*15 (01, 04-07, 12, 14, 19, 20, 24, 25, 26N, 27, 32, 33, 
34 and 35); group II: HLA-B*15 (02, 08, 11, 15, 28 and 
30); group III: HLA-B* 1503/4802; group IV: HLA- 
B*1521. Groups I and II occurred frequently (allele 
frequency = 8.0 and 2.5%), and thus we optimized the 
polymerase chain reaction-single-stranded conformation 
polymorphism (PCR-SSCP) method to identify HLA- 
B*15 subtypes of groups I and II. Eighty samples of 
DNA carrying HLA-B* 15 from 300 healthy unrelated 
individuals were tested. B*1502 (52.5%) and B*1525 
(13.8%) were the most common subtypes found in NET. 
They also showed strong linkage disequilibrium with 
HLA-Cw and heterogeneity of HLA-A, DR, DQ haplo- 
types. Although limited conclusions can be drawn from 
this study because of the small number of DNA references 
used, the baseline data will be useful in the selection of 
common HLA-B* 15 alleles when subtyping for unrelated 
donor transplantations. 

Introduction 

Human leukocyte antigen (HLA) B15 represents a large, 
genetically diverse group of HLA-B alleles. It was first 
described and its heterogeneity recognized at the Fourth 
International Histocompatibility Workshop in 1970 (4th 
IHW) (Alben et aL, 1970; Thorsby &; Kissmeyer-Nielsen, 
1970; Thorsby etaL, 1970). At the 9th IHW, the heter- 
ogeneity of HLA-B15 was summarized and discussed by 

" Department of Clinical Immunology, Faculty of Associated Medical 
Sciences, t Center for Research and Development of Medical Diagnostic 
Laboratories, t Faculty of Associated Medical Sciences, and § Blood 
Transfusion Center, Faculty of Medicine, Khon Kaen University, Khon 
Kaen, Thailand. 

Received 26 June 2002; revised 27 August 2002; 
accepted 1 1 November 2002 

Correspondence: Amornrat V. Romphruk, Blood Transfusion Center. 
Faculty of Medicine, Khon Kaen University, Khon Kaen 40002, Thailand. 
TetyFax: 6M3-246942; E-mail: arunrant@lcku.ac.th 



Chandanayingyong et al. (1984) and Cambon-Thomsen 
et al, (1984). At the 10th IHW, seven serological subtypes 
of B15 were identified and five were recognized by the 
World Health Organization (WHO) Nomenclature 
Committee as distinct specificities: B62, B63, B75, B76 
and B77 (Albert et al., 1989). Many new HLA-B15 
related alleles have since been discovered in different ethnic 
groups throughout the world, especially among South- 
east Asians (Alonso et a/,, 1983; Cambon-Thomsen et ai, 
1984; Albert etaL, 1989; Parham etal., 1989; Pohla 
etai, 1989; Nisperos etaL, 1991; Choo etaL, 1991, 
1993; Little & Parham, 1991; Belich etaL, 1992; San- 
tamaria et aL, 1993; Hildebrand et aL, 1994; Wang et aL, 
1997). The correlations of HLA-B15 molecular structures 
with serological specificities and antigenic properties have 
been discussed (Pohla etaL, 1989; Hildebrand etaL, 
1994), but the serological assignments of B*15 alleles to 
B62, B63, B75, B76 and B77 are complicated by cross- 
reactivities among the antigens encoded by these alleles 
(Hildebrand et aL, 1994; Lin et aL, 1996; Sterner etaL, 1997; 
Marsh et aL, 2001 ). Although the heterogeneity of HLA-B15 
antigens has been described in both serological and molecular 
analyses, the serological complexities of this group and 
the corresponding molecular structures are still unclear. 

In the Thai population, HLA-B15 antigens have been 
investigated in distinct ethnic groups such as the popula- 
tion of central Thailand (Bangkok) (Chandanayingyong 
etaL, 1997), Dai Lue (Chandanayingyong etaL, 1997), 
the Thai and Thai-Chinese population (Imanishi et aL, 
1991), the population of northern Thailand (Fongsatikul 
etaL, 1997) and that of southern Thailand (Chiewsilp 
etaL, 1997). The most common serologically defined 
antigens were B75 and B62. 

We have studied the distribution of HLA-B* 15 alleles in 
300 samples from the population of north-eastern Thai- 
land (NET) based on the 12th IHW PCR-amplification 
refractory mutation system (ARMS). Six groups of 
B*15 could be defined by this protocol. Four groups of 
B*15 were identified in this population: B*15 (01, 04-07, 
12, 14, 19, 20, 24, 25, 26N, 27 and 32-35) [allele fre- 
quency (AF) = 3.67%], B*15 (02, 08, 11, 15, 28 and 30) 
(AF= 8.67%), B*1503/4802 (AF = 0.5%), and B*1521 
(AF = 0.5%) (Romphruk et aL, manuscript in preparation). 
Each group was carried by different HLA class I and class 
II haplotypes. In spite of these differences, the biological 
significance of the genetic heterogeneity of HLA-B15 
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requires to be elucidated. Before this information can be 

collated, the subtypes of HLA-B15 need to be identified. 
Thus, our aim was to identify HLA-B*15 subtypes in 
NET by PCR-single-stranded conformation polymor- 
phism (PCR-SSCP). 

Materials and methods 

Samples 

A total of 80 HLA-B15 positive samples were obtained from 
300 unrelated healthy individuals from north-eastern 
Thailand (Romphruk etal, manuscript in prepara- 
tion). All of these samples were typed for HLA-B*15 
alleles using the 12th IHW PCR-ARMS protocol. Four 
groups were differentiated, but only two were studied by 
PCR-SSCP, namely: group I (« = 22), B*15 (01, 04-07, 
12, 14, 19, 20, 24, 25, 26N, 27 and 32-35) and group II 
(n = 52), B*15 (02, 08, 11, 15, 28 and 30). The DNA 
references for B*15 aUeles (namely B*1501, 1502, 1503, 
1505, 1507, 1508, 1510, 1511, 1512, 1513, 1516, 1518, 
1520, 1524, 1525, 1527, 1528 and 1532) were provided 
by the Japanese Red Cross Central Blood Center, Tokyo, 
Japan and the Australian Red Cross Blood Service, 
Melbourne, Australia. 

Primers 

A two-step PCR was performed (Fig. 1) using the same 
nucleotide primer sequence as Lin et al (1996). In the first 
step, the HLA-B*15 gene fragment was amplified from 
genomic DNA by PCR using the group-specific primers 
BEX2-1 and BCT. The amplified fragment was about 
742 bp in length encompassing exon 2, intron 2 and most 
of exon 3, A second PCR was then performed separately 
to amplify exon 2 (215 bp) and exon 3 (207 bp) using primer 
pairs BEX2-1/BEX2-FC and BEX3-1/BCT, respectively. 



PCR amplifications 

PCR amplifications were performed in a 10-p.l reaction 
mixture. In the first PCR, the mixture consisted of 100 ng 
of genomic DNA as template, 0,5 |im of each primer, 
200 ^IM of each dNTP, 50 mM Tris-HCl (pH 8,8), 10 mM 
(NH4)2S04, 1.0 mM MgCl2 and 0.5 units of Taq DNA 
polymerase (Promega, Madison, WI or MHC Research 
Group, Khon Kaen University (Leelayuwat et al.y 1997)). 
After preheating at 94 ®C for 2 min, 35 PCR cycles were 
performed, followed by final extension at 72 °C for 
10 min. Each amplification cycle consisted of 1 min each 
of denaturing (94 °C), anneahng (54 °C) and extension 
(72 °C). The PCR was carried out in a Gene Amp PCR 
system 9600 (Perkin Elmer Corporation, Norwalk, CT). 

The reaction mixture for the second PCR was the same 
as that for the first PCR, except that we used 1 fil of the 
first PCR product as the template, 0.8 mM MgClj for 
exon 2 and 1.2 mM MgCli for exon 3. The ampHfication 
was performed for 20 cycles, after preheating to 94 °C 
for 2 min and final extension of 10 min at 72 **C. Each 
amplification cycle consisted of denaturing (94 °C) for 30 s, 
annealing (58 °C) for 30 s, and extension (72 °C) for 1 min. 
The presence or absence of PCR product was determined 
by gel electrophoresis. Five microlitres of PCR product 
was run on a 1 % agarose gel containing 0.5 \Lg ml"' of 
ethidium bromide and visualized under UV illumination. 

SSCP fragment analysis 

The exon 2 and 3 fragments ampUfied by the second 
PCR were analysed separately using the SSCP method 
(Lin etaL, 1996; Maruya et ai, 1996). Elearophoresis 
was carried out in a Miniprotein II polyacrylamide gel 
(Mini-PROTEAN II, Biorad, Hercules, CA) (Table 1). The 
single-stranded DNA fragments separated in the gel were 
detected by silver staining. 
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Rgure 1. The strategy of SSCP fragment 
analysis of axons 2 and 3 in HLA-8*1 5 alleles. 
Prinners BEX2-1 and BCT are group-specific 
primers. BEX2-1/BEX2-FC and BEX3-1/BCT 
amplified the exon 2 and exon 3 seperately. 
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Table 1. SSCP conditions for groups I and II of the B*1 5 alleles 
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39: 1 
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exon 3 
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50: 1 


19. 20.25. 26N. 27. 32-35) 














Group II 


exon 2 
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150 


150 


10% 


50: 1 


BM5(02. 08, 11. 15. 28. 30) 


exon 3 


4 


150 


60 


7% 


50:1 




Rgure 2, PCR-SSCP analysis of HLA-BI 5 
group 1(1501, 1505. 1507, 1512, 1520. 1524, 
1525, 1527. 1532). Lanes 1 (8*1525), 2 
(B*1532),3(B*1507),4(B*1520), 5(B* 1512). 
6 IB' 1524), 7 (B*1527). 8 (B* 1501) and 
9 (B* 1505) were distinguished using exon 3 
fragment analysis. Lanes 10 (B*1501) and 1 1 
{BM524)byexon 2. 
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RgureS. PCR-SSCPanalysisof HLA-BI 5 Group I! (1502. 1508. 1511, 
1528). Lanes 1 (B*1502) and 2 (8*1508/151 1/1528) were done using 
exon 3 fragments analysis. Lanes 3 (B*1508), 4(B*151 1) and 
5(B*1528) by exon 2. 



Reproducibility of the test 

To ensure reproducibility, the DNA references were tested 
within and between runs for both PGR amplification and 
SSCP fragment analysis. 

Results 

SSCP fragment analysis 

Most of alleles in group I had nucleotide sequence differ- 
ences in exon 3, so analysis was used to discriminate alle- 
les B*1505, 07, 12, 20, 25, 27 and 32. The other two 
alleles, B*1501 and B*1524, were defined by the exon 2 
SSCP analysis (Fig. 2). In group II, exon 3 SSCP analysis 
discriminated B*1502 from B*1508, 1511 and 1528 and 
exon 2 discriminated B*1508, B*1511 and B*1528 (Fig. 3). 

Reproducibility of the test 

The SSCP analysis of DNA references for exon 2 for both 
groups I and II showed identical banding patterns Vithin' 



and 'between' runs. The exon 3 analysis showed slight dif- 
ferences in banding patterns on the 'between run' assay, 
but these alleles were easily defined. Thus, the DNA 
references were included in every SSCP analysis. 

B*15 alleles in NET 

Eight different patterns of HLA-B* 15 alleles (groups I and 
II) were detected in NET (Table 2). The common alleles 
were B*1502 (52.5%), B*1525 (13.8%) and B*1528 
(5.0%). Three samples of group I [B*15 (01, 07)?] could 
not be clearly defined because extra bands froni the other 
HLA-B allele in these samples could not be excluded. One 
of the 22 samples from group I and two of the 52 in group 
II were undefined because they gave different banding 
patterns from the DNA references. They might have been 
B*15 (04, 06, 14, 19 or 33-35) or BM5 (15 or 30) for 
which we have no DNA references. Five samples from 
group I and three from group II could not be analysed 
because DNA samples were not available. The com- 
parison of B*15 alleles in NET and other populations is 
presented in Table 3. 

Haplotype analysis 

Table 4 summarizes the possible association of HLA- 
B*15 alleles with other HLA alleles. The data for other 
loci are from a previous study (Romphruk et aiy manu- 
script in preparation). Only four haplotypes in Table 4 
could be confirmed by family segregation analysis. Most 
HLA-B*1502 haplotypes (92.85%) carried HLA-Cw*08 
but carried diverse alleles at HLA-A, -DR and -DQ. All 
HLA-B*1525 and B*1521 haplotypes carried HLA- 
Cw*04 but carried diverse alleles at HLA-A, -DR and 
-DQ. In addition, associations of B*1528 and B*1532 
with Cw*07 (01-03) were observed. 
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B*15* 


B* 15 alleles 


n (%) in = 80) 


Group 1 


1501 


— 




1505 


— 




1507 






1512 


— 




1520 






1524 


— 




1525 


1 1 (13.8%) 




1527 


— 




1532 


2 (2.5%) 




15(01,07)? 


3 (3.7%) 




15 (04, 06. 14, 19, 33, 34, 35) 


1 (1.3%) 




Not done 


5 (6.2%) 


Group II 


1502 


42 (52.5%) 




1508 






1511 


1 (1.3%) 




1528 


4 (5.0%) 




15(15, 30)? 


2 (2.5%) 




Not done 


3 (3.7%) 


Group III 


1503/4802 


3 (3.7%) 


Group IV 


1521 


3 (3.7%) 


Group V 


15(13, 16. 17) 


0 


Group VI 


15(09. 10, 18, 23) 


0 



Table 2. 8*15 alleles in HLA-B1 5 positive 
samples in = 80) from 300 healthy individuals 
from north-eastern Thailand 



' The 815 groups were defined according to the 12th IHW PCR-ARMS. 



Table 3. Percentage of HLA-B* 15 subtypes in various populations 



B*15 


NET (n = 80) 


Vietnamese* 
(n = 43) 


Caucasian* 
in = 70) 


Ni-Nanuatu" Bubi* 
(n = 37) (n = 30) 


Koreans' 
in = 237) 


Japanese' 
(n = 53) 


1501 


3.75' 


20.9 


85.7" 


5.4 


35.0" 


30.2" 


1502 


52.5 


55.8 






1.5" 


0.9" 


1503 


3.7" 






53.3" 




ND 


1506 


ND 






89.19" 




ND 


1507 






1.4 




4.6 


4.7 


1509 






1.4 






ND 


1510 








40.0" 




ND 


1511 


1.3 








7.6" 


2.8 


1512 




14.0" 








ND 


1515 


2.5* 








ND 


ND 


1516 








10.0 




ND 


1517 






5.7 




0.002 


ND 


1518 


ND 




7.1 




3.0 


10.4 


1521 


3.7 


2.3 




8.1 




ND 


1525 


25.6 


14.0 






0.002" 


ND 


1527 




2.3 






1.5 


0.9 


1528 


5.0 










ND 


1532 


2.5 








ND 


ND 


1538 


ND 


ND 


ND 


ND ND 


0.004 


ND 



Blank = not found. ND = not done/no data. ' Bl 5 (01 or 07); " B1 503/4802; B1 5 (1 5 or 30); " significant difference compared to NET (corrected P < 0.05); 
• data from Bernardo etal. (1998); ' data from Lee etal. (2000); « data from Saito ef al. (2000). 



Discussion 

HLA-B* 15 is a diverse group of alleles of the HLA-B locus 
and is found throughout the world, suggesting that its 
lineage is as old as Homo sapiens (Hildebrand^t ai, 1994). 
By DNA sequencing of exons 2 and 3, 68 B*15 alleles 
were defined (Marsh et al., 2001). B*15 alleles associated 
with specific ethnic groups have been described, including 



B*1518 in Asian Indians and B*1521 in Australian 
Aborigines (Lienert et aL, 1995). 

We used PCR-SSCP to define the B*15 allele in NET 
This technique depends on the conformation of single- 
stranded DNA to produce the different banding patterns, 
so DNA references are needed for comparisons. However, 
the references themselves varied under different conditions. 
Unfortunately, we did not have all the DNA references. 
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Table 4. Possible associations of B*15 alleles with other HLA loci 



s* 


A 


Cw* 


DRB1* 


DQBl* 


% 


1502 




08 






92.85 (39/42) 




n 


08 


1202 


0301 


26.19(11/42) 




24 


08 


1501 


0601 


26.19(11/42) 




n 


08 


1502 


0501 


4.76 (2/42) 


1525 




04 






100(11/11) 






04 


1106 


0301 


27.27 (3/11) 






04 


1502 


0501 


18.18(2/11) 






04 


1502 


0502 


18.18(2/111 


1521 


34/66 


04 






100 (3/3) 


1528 




07(01-03) 






50 (2/4) 


1532 




07(01-03) 






100(2/2) 



SO only nine of the 17 alleles for group I [B*15 (01, 05, 

07, 12, 20, 24, 25, 27 and 32)] and four of the six for 
group II [B*15 (02, 08, 11 and 28)] were included in our 
study. Specific primers can amplify other alleles at the B 
locus, such as B*4601, B*5801, B*38, B*39 and B*35 
(Lin etaL, 1996), so, in addition to the two alleles per 
locus in individual samples, the SSCP pattern revealed 
other bands which were not from B*15 alleles. Thus, 
including other known alleles present in the sample in 
the SSCP fragment analysis will help to clarify the B* 15 
bands. 

HLA-B 15 allelic types and their frequencies have 
previously been shown to vary in different populations 
(Marsh et ai, 2001; Lee et ai, 2000). Frequencies of 
some of the most common B*15 alleles have been 
determined: B*1502 in NET (52.5%) and Vietnamese 
(55.8%); B*1501 in Caucasians (85.7%), Koreans (35.0%), 
and Japanese (30.2%); B1506 in Ni-Vanuatu (89.2%); 
B*1503 and B*1510 in Bubi (53.3 and 40.0%). The dis- 
tribution of the B*15 alleles in NET is similar to that in 
the Vietnamese, except for B* 1512. Compared to other Thai 
groups (Thais and Thais-Chinese, Thais from Central Thai- 
land and Dai Lue), B*75 (encoded by B*1502, 1508, 
1511 and 1521), B62 (encoded by B*1501, 1504-1507, 
1515, 1520, 1524, 1525, 1527, 1530, 1532, 1545, and 
1548) and B77 (encoded by B*1513) were found more 
frequently (Chandanayingyong etal.^ 1997; Imanishi 
etai, 1991). It is possible that B75 was encoded by 
B*1502 and that most of B62 was encoded by B*1525 in 
NET. In contrast, B*1501 was the predominant allele for 
the B62 antigen in Koreans (166 in 198 samples) but 
B*1525 was rare (1 in 198 samples) (Lee etaL, 2000). 
B*1513 was not found in the present study, indicating 
that this allele is rare in north-eastern Thailand. In a sep- 
arate study, we found this allele in only one in 240 cases 
of cholangiocarcinoma in patients from north-eastern 
Thailand (Romphruk et ai, manuscript in preparation). 
This confirms our observation that the distribution of 
MHC alleles in the population of north-eastern Thailand 
is different from that in the population of Central Thailand, 
which shows a high degree of Thai-Chinese admixture 
(Romphruk et aL, 1999), 

HLA-B15, combined with B4601, is the most common 
(22-30%) HLA-B allele in south-east Asians, including 



the Chinese. Hildebrand et al (1994) and Lin et al (1996) 
proposed a phylogenetic tree for 26 HLA-B 15 alleles and 
a network of structure relationships among these alleles. 
B4601 may have evolved from B*1501 — the putative 
original B*15 allele (Zemmour etal, 1992; Hildebrand 
et al, 1994; Lin et al, 1996;), B*4601 and B*1501 differ 
in the short codon position 66-76, which may be the 
result of segmental conversion between B*1501 and 
Cw*01. Thus, HLA-B*4601 in Asians may have origi- 
nated from B*1501 via gene conversion. This hypothesis, 
combined with selective pressure, may explain the low 
frequency of B* 1501 and the high frequency of B*4601 
in these populations. 

The possible associations of B*15 with other loci are 
presented in Table 4. The possible B* 15-Cw* association 
showed several different patterns. B*1502 was strongly 
associated with Cw*08, B*1525 with Cw*04, and 
B*1528 and B*1532 with Cw*07 (01-03). B*1502 and 
B*1525 were more diverse in their associations with 
HLA-A, Cw, DR and DQ than other B*15 alleles. This 
suggests that these two alleles may be relatively old alleles 
in the Thai population, similar to B*1501 and B*1518 in 
Koreans. The diversity of haplotypes is likely to be the 
result of recombination between HLA-B and other loci 
during the long period of haplotype formation in the pop- 
ulation (Lee et al, 2000). The strong linkage disequilib- 
rium between B*15 alleles and HLA-Cw will be useful in 
discriminating the B15 subgroups, especially where B*15 
subtyping is not available. 

In conclusion, we defined B*15 alleles and the hetero- 
geneity of haplotypes carrying HLA-B15 in NET. High- 
resolution typing will be necessary to select unrelated 
donors for organ and bone marrow transplantation. The 
data obtained in this and similar studies will facilitate the 
selection of common B*15 alleles for subtyping in these 
populations. 
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Abstract: MICA or PERBll.l is a polymorphic major histocompatibility 
complex (MHC) class I-related gene located 46 kb centromeric of the HI.A-B 
gene in the HLA class I region. It is expressed mainly in gut epithelial cells, 
keratinocytes, endothelial cells, fibroblasts and monoc>tes, and is upregu- 
lated by heat stress. MICA has been found to interact with yb T cells, ap 
CD8 ' and natural killer (NK) cells bearing the NKG2D/DAP10 receptor. The 
MICA gene displays a high degree of polymorphism with at least 54 alleles. 
In the present study, polymoi-phic exons 2, 3 and 4 of the MICA gene were 
analyzed using sequencing based typing (SBT) in 255 unrelated healthy 
northeastern Thais. Thirteen previously reported MICA alleles were de- 
tected M]CA*008, *010, *002 and *019 were highly predominant with the 
allele frequencies of 21.4%. 18.2%, 17.6% and 15.3%, respectively Five of 
these 13 MICA alleles show significantly different frequendes from those 
of the Japanese and Caucasian populations. Interestingly, M1CA052, which is 
a very rare allele in other populations, was prevalent with the allele fre- 
quency of 8.2%, mainly on the HLA haplotype cai-iying HLA-B*13 in this 
population Strong linkage disequilibria were observed between MICA and 
HLA-B, as similarly observed in other populations, namely MICA*010- 
B*4601. MICA052-B*13, MICA*002.B*5801. and MICA*019-B*15 (1502. 
1508. 1511, 1515, 1528. 1530). A large variety of three-locus (MICA - HLA- 
B - ra.A-Cw) and six-locus aH^A-DQBl - HLA-DRBl - MICA - ULA-B - 
HLA-Cw - HLA-A) haplotypes were recognized in the northeastern Thai 
population. This is the first report on MICA allelic distribution in Southeast 
Asian populations. These data will provide the important basis for future 
analyses on the potential role of the MICA gene in disease susceptibility 
and transplantation matching in Southeast Asian populations. 

The telomeric region of the human major histocompatibility complex 
(MHC) includes the MHC class I chain-related (MIC) (1) or PERBll (2) 
gene family consisting of seven members: MICA (PERBll.l), MICB 
(PERB11.2), Mice (PERB1L3), MICD (PERB1L4), MICE (PERB11.5), 
MICF and MICG (3, 4). Of these, only MICA and MICE are expressed, 
encoding cell surface glycoproteins of 383 amino acids. The MICA 
gene is located approximately 46 kb centromeric of HLA-B. MICA 
molecules are composed of three extracellular domains (al, ot2 and 
otS), a transmembrane segment and a carboxy-terminal cytoplasmic 
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tail similar to other HLA class I molecules. However, MICA shares 

only an average of 21, 19 and 34% amino acid identity in the al, a2 
and (x3 extracellular domains, respectively, with classical HLA class 
I proteins. Further, MICA and HLA class I gene products differ re- 
garding tissue distribution, regulation of expression, association with 
P2-microglobuIin and immunological function. 

MICA has a divei'se tissue distribution. It is transcribed in fibro- 
blast, epithelial cell lines, gastrointestinal epithelium, freshly isolated 
keratinocytes, endothelial cells and monocytes but not in B lympho- 
cyte (5, 6). MICA is also expressed in some lung, breast, kidney, ovary, 
prostate and colon carcinoma (7). Expression of the gene is not 
affected by y-interferon stimulation but is regulated by promoter 
heat-shock elements and strongly induced upon cell stress, thus 
MICA may function as an indicator of cell stress (8, 9). MICA is recog- 
nized by a subset of y5T cells that predominate in intestinal and other 
epithelial sites (8, 9), which suggests that MICA facilitates y5T-cell 
detection of infected, damaged, or otherwise stressed epithelial cells, 
thus forming part of the innate immune system (10). Recently, MCA 
has been reported to be a ligand for an activating receptor, >^KG2D, 
which is expressed on most natural killer (NK) cells, CD8+aPT cells 
and yST cells (11). This interaction might activate NK-cell and T-cell 
responses against MICA-bearing tumors. 

The MICA gene displays a high degree of genetic polymoi-phism 
in exons 2, 3 and 4, amounting to at least 54 alleles (12-14). In 
contrast to classical class I genes, the polymorphic sites are distrib- 
uted throughout the domains but outside the peptide-binding 
groove of the MHC class I molecules (12). They are defined by a 
total of 40 nucleotide substitutions, 30 of which are nonsynony- 
mous; 5 of 8 in exon 2, 15 of 16 in exon 3 and 10 of 16 in exon 4 
(13). Genetic polymorphism was also recognized in the transmem- 
brane region including indels and variable number of the GCT trip- 
let repeats (4). However, information on the MICA diversity in hu- 
man ethnics is limited, being derived mainly fi-om the Caucasian 
and Japanese populations. In the present study, we have used se- 
quencing based typing (SBT) to define MICA alleles in the native 
Northeastern Thai (NET) population. The results were compared 
with those from the Caucasian and Japanese populations. Linkage 
disequilibria among MICA and HLA alleles were also investigated 
to identify representative HLA haplotypes in the NET population. 



Material and methods 

Genomic DNA 

Peripheral blood cells were collected from 255 unrelated healthy 
Northeastern Thai individuals. All of them were interviewed for 



their ancestors. Their families are living in the northeast of Thai- 
land for at least two generations. Genomic DNAs were extracted 
from buffy coat by the salting-out method (15). Their HLA class I 
(HLA-A, -B and -C) and class II alleles (HLA-DRBl and -DQBl) 
were typed by the polymerase chain reaction-sequence-specific 
primer (PCR-SSP) technique (Romphruk et al., in preparation). Gen- 
omic DNA from an HLA homozygous tissue culture cell line 
(HTLC), BM15 with MICA*004 was employed as a standard DNA 
in the SBT protocol. 

PCR amplification and sequence determination of the MICA 
gene 

A 2.2-kb MICA gene fragment from exon 2 to exon 5 was amplified 
by a pair of PCR primers. 5MICA and 3MICA (12). PCR conditions 
using these primers were desanbed previously (16) with some modi- 
fications. Namely, PCR was carried out in a 50-|il reaction using 
100-200 ng of DNA and 0.5 units AmpliTaq™ DNA polymerase in 
the GeneAmp® PCR System 9700 (Applied Biosystems, Foster City, 
CA, USA). The cycling conditions were 95°C for 15 min, followed 
by 35 cycles of 95°C 30 s. 6rC for 1 min, 72°C for 2 min and final 
elongation at 72''C for 10 min. 

Cycle sequencing 

The PCR products were purified using the QIAquick PCR Purifi- 
cation Kit (Qiagen, Germany) and subjected to sequence determi- 
nation using the ABI PRISIVP*^ BigDye Terminator Cycle Sequenc- 
ing Kit with AmpliTaq DNA polymerase, FS (Applied Biosystems). 
Each of exons 2, 3 and 4 was separately sequenced on both DNA 
strands using six sequencing primers as described by Katsuyama 
et al. (16). The products were purified by Centri-Sep'^'^ spin columns 
(Princeton Separations, Adelphia, NJ, USA). The reactions were run 
on the ABI 377 or ABI 310 Sequencing System (Applied Biosys- 
tems) and the results were analyzed using the ABI Sequence Analy- 
sis and Navigator Softwares, Raw sequencing data were manually 
inspected for confirmation and alleles were assigned by alignment 
with known 54 MICA alleles of the extracellular domains (13, 14). 

Statistical analysis 

Allele frequency (AF) was calculated using the following formula: 
%AF=sum of each individual allele/2N xlOO. Haplotype frequency 
(HF) and linkage disequilibrium (D) were calculated by the maxi- 
mum likelihood method assuming the Hardy- Weinberg law. This 
value indicates an association between the two loci. Significance 
of two-locus association was assessed by the Chi-square test. The 
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significance of the distribution of alleles between NET and other 
populations was tested by Pc value (corrected P-value) after the Chi- 
square method. 

Results 

MICA allele frequencies in Northeastern Thais 

The distribution of MICA alleles was determined by sequencing 
based typing in 255 unrelated healthy Northeastern Thai individ- 
uals (NET). Thirteen alleles were recognized in this population. The 
allele fi*equencies (AF) of the MICA gene are listed in Table 1, The 
common alleles were MICA*008 (21.4%), MICA*010 (18.2%), 
MICA*002 (17.6%) and 1VGCA*019 (15.3%). In contrast, MICA*020 
(0.4%) and *016 (0.2%) were infrequent. Most of the *^lank" allele 
cases revealing the frequency of 7.8% possibly represent a MICA 
homozygote, but some or a significant fraction of them may corre- 
spond to a MICA deletion allele linked to HLA-B*4801 in Asian 
populations (17, 18). In fact, B*4801 was not so frequent, but was 
recognized with the allele frequency of around 0.5% in NET 
(Romphruk et al., unpublished data). 

Comparison of MICA allele frequencies among NET, 
Japanese and Caucasians 

The allelic distribution of the MICA gene in NET was compared 
with those of the Japanese (16) and Caucasian (19) populations 

Allele frequencies of MICA In 
255 Northeastern Thais (NET) 



MICA 


n 


% AF 


002 


90 


17.6 


004 


18 


3.5 


008 


109 


21.4 


009 


12 


2.4 


010 


93 


18.2 


012 


16 


3.1 


016 


1 


0.2 


017 


11 


2.2 


018 


31 


6.1 


019 


78 


15.3 


020 


2 


0.4 


026 


7 


1.4 


052 


42 


8.2 



Table 1 
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Comparison of MICA alleles of NET with those of Japanese and Caucasians 





%AF 








NET 


Japanese (16)** 


Caucasians (19)** 


MICA 


n--255 


n = 130 


n-242 


001 


0.0 


0.0 


3.0* 


002 


17.6 


14.6 


13.0 


004 


3.5 


9.2* 


6.0 


007 


0.0 


1.2 


2.0* 


008 


21.4 


30.8* 


55.0* 


009 


2.4 


16.5* 


3.0 


010 


18.2 


10.8* 


5.0* 


Oil 


0.0 


0.0 


2.0 


012 


3.1 


12.3 


2.0 


016 


0.2 


0.0 


2.0 


017 


2.2 


0.0 


4.0 


018 


6.1 


0.0* 


2.0* 


019 


15.3 


3.5* 


0.0* 


020 


0.4 • 


0.0 


1.0 


026 


1.4 


0.0 


0.0 


051 


0.0 


0.0 


2.0* 


052 


8.2 


0.0* 


0.0* 


* significantly different (Pc<0.05) 






•* reference no. 








Table 2 



(Table 2). As a result, it must be most notable that MICA052 
(MICA052 has not be officially recognized and so * is not attached 
between MICA and 052) (14) was considerably prevalent with the 
frequency of 8.2% in NET, but absent in Japanese or (^ucasians. 
Further, the significantly high frequencies of MICA*010 (the second 
most frequent allele in NET), MICA*018 and MICA*019 were found 
in NET as compared to the Japanese and Caucasian populations by 
the corrected P-value test. On the other hand, the frequencies of 
MICA*008 were significantly decreased in NET as compared to 
those of both populations, however, MICA*(X)8 was the most fre- 
quent allele amongst all the three populations including NET. Tn 
addition, the significantly low frequency of MICA*004 and 
MICA*009 was observed in NET, only when compared to the Ja- 
panese population. There was no MICA*001, MICA*(X)7 or 
MICA051 (MICA051 has not be officially recognized) (14) allele in 
NET, which was present with the low frequency in the Japanese 
and Caucasian populations. 

Linkage disequilibrium and liapiotypes between iVIiCA and 
HU B in NET 

Linkage disequilibrium parameters (delta values) and haplotype fre- 
quencies (HF) between MICA and HLA-B alleles were calculated in 
NET. Strong associations were observed between these two loci as 
previously reported in the Japanese population (16). MICA - HI-A- 
B associations with statistical significance at the P level of less than 
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Significant MICA and HLA-B associations In 255 NET 



MICA HLA-B 



002 


35 


20 


3.9 


2.3 


21.5 


002 


38 


14 


2.8 


1.7 


17.5 


002 


5801 


49 


9.6 


7.1 


94.8 


004 


44 


17 


3.3 


3.1 


224.8 


008 


27 


26 


5.1 


3.3 


27.5 


008 


4001 


33 


6.5 


4.6 


45.5 


008 


07 (02-07) 


13 


2.5 


1.1 


4.7 


008 


39/6701 


21 


4.1 


2.3 


14.4 


008 


40 (02,04-06) 


16 


3.1 


1.8 


16.9 


009 


51/5201 


10 


2.0 


1.9 


99.3 


010 


4601 


75 


14.7 


11.6 


160.6 


010 


15a 


18 


3.5 


2.4 


17.0 


012 


5401 


4 


0.8 


0.8 


60.7 


012 


55/56 


11 


2.2 


2.0 


99.0 


017 


57 


11 


2.2 


2.1 


232.8 


018 


18 


30 


5.9 


5.1 


158.3 


019 


1521 


8 


1.6 


1.2 


16.3 


019 


15a 


13 


2.5 


1.4 


6.9 


019 


15b 


46 


9.0 


6.5 


87.5 


026 


07 (02-07) 


6 


1.2 


0.9 


52.0 


052 


13 


41 


8.0 


7.6 


215.4 



15a- 15 (01, 04-07. 12, 14. 19. 20. 24. 25. 26N. 27. 32-35) 
15b«15 (02. 08. 11. 15. 28. 30) 



Tables 



0.05 corresponding to a delta- value of equal to or more than 0.8 are 
listed in Table 3 with 21 representative haplotypes. The most com- 
mon haplotype was MICA*010-B*4601 (14.7%), followed by 
MICA*002-B*5801 (9.6%), MICA*019-B*15(02,08.11,15,28,30) (9,0%) 
and MICA052-B*13 (8.0%). The NET population was characterized 
by a large variety of associations of MICA with HI.A-B, such as 
MICA*002, MICA*008, MICA*010. ]VIICA*012 and MICA*019 
showing linkage disequilibria with multiple HLA-B alleles as shown 
in Table 3. 

MultMocus haplotypes in NET 

Frequencies of three-locus (MICA - HLA-B - HLA-Cw) and six- 
locus (HLA-DQBl - HLA-DRBl - MICA - H1.A-B - HLA-Cw - 
HLA-A) haplotypes were calculated by incorporating the above 
MICA - HLAB haplotypes into the known HLA class I and class n 
haplotypes in NET, and those with more than 0.4% (more than two 
individuals carrying the same haplotype) are listed in Tables 4 and 
5, respectively. They revealed the presence of diversed and con- 



Posslble MICA - B - Cw haplotypes in 255 NET 



008 


07 (02-07) 


07 (01-03) 


11 


2.2 


026 


07 (02-07) 


15 (02.03.05) 


5 


1.0 


052 


13 


0304 


22 


4.3 


052 


13 


04 


18 


3.5 


008 


13 


0602 


5 


1.0 


019 


15a 


04 


9 


1.8 


010 


15a 


07 (01-03) 


6 


1.2 


019 


15b 


04 


4 


0.8 


019 


15b 


08 


36 


7.1 


010 


15b 


07 (01-03) 


2 


0.4 


019 


1521 


04 


7 


1.4 


008 


18 


07 (01-03) 


4 


0.8 


019 


18 


07 (01-03) 


5 


1.0 


018 


18 


0704 


28 


5.5 


018 


18 


04 


2 


0.4 


008 


27 


0304 


17 


3.3 


008 


27 


07 (01-03) 


2 


0.4 


019 


27 


1202 


3 


0.6 


002 


35 


04 


13 


2.6 


002 


38 


07 {01-02) 


10 


2.0 


002 


39/6701 


07 (01-03) 


2 


0.4 


008 


39/6701 


07 (01-03) 


18 


3.5 


008 


40 (02,04^6) 


15 (02.03.05) 


5 


1.0 


008 


40 (02,04-06) 


0304 


2 


0.4 


008 


4001 


0303 


3 


0.6 


008 


4001 


0304 


7 


1.4 


008 


4001 


04 


6 


1.2 


008 


4001 


07 (01-03) 


14 


2.8 


004 


44 


07 (01-03) 


16 


3.1 


010 


4601 


01 


73 


14.3 


020 


5001 


0602 


2 


0.4 


019 


51/5201 


07 (01-03) 


2 


0.4 


002 


51/5201 


14 


6 


1.2 


009 


51/5201 


14 


5 


1.0 


009 


51/5201 


15 (02.03.05) 


3 


0.6 


012 


5401 


01 


4 


0.8 


012 


55/56 


01 


6 


1.2 


012 


55/56 


1202 


2 


0.4 


002 


55/56 


1203 


2 


0.4 


017 


57 


0602 


11 


2.2 


002 


5801 


0302 


42 


8.2 



I5a»l5 (01. 04-07, 12. 14. 19. 20, 24. 25. 26N. 27. 32-35) 
I5b=l5 (02, 08. 11. 15. 28. 30) 

Table 4 



n %Hf Delta Chi-square MICA HLA-B HLA-Cw n %HF 
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Possible MICA and 5- locus HLA haplotypes in 255 NET Table 5 



DQBl 


DRBl. 


MICA 


B 


Cw 


A 


n 


%HF 


0501 


1001 


026 


07 (02-07) 


15 (02,03,05) 


29 


4 


0.8 


0501 


1502 


052 


13 


04 


11 


4 


0.8 


0502 


1602 


052 


13 


0304 


11 


5 


1.0 


02 


07 


008 


13 


0602 


30 


5 


1.0 


0502 


1502 


010 


15a 


07 (01-03) 


34/66 


4 


0.8 


0301 


1202 


019 


15b 


08 


11 


10 


2.0 


0601 


1501 


019 


15b 


08 


24 


6 


1.2 


0601 


1502 


019 


1521 


04 


34/66 


6 


1.2 


0501 


1502 


019 


18 


07 (01-03) 


0203 


3 


0.6 


0501 


1502 


018 


18 


0704 


0203 


5 


1.0 


0501 


1502 


018 


18 


0704 


24 


8 


1.6 


0301 


1202 


018 


18 


0704 


24 


5 


1.0 


0301 


1101 


008 


27 


0304 


24 


2 


0.4 


0301 


1202 


008 


27 


0304 


24 


3 


0.6 


0501 


1502 


008 


27 


0304 


24 


6 


1.2 


0501 


1502 


002 


38 


07(01-03) 


11 


2 


0.4 


0501 


1502 


002 


38 


07 (01-03) 


24 


2 


0.4 


0302 


0405 


008 


39/6701 


07 (01-03) 


0203 


5 


1.0 


02 


07 


004 


44 


07 (01-03) 


33 


13 


2.5 


0303 


09 


010 


4601 


01 


11 


6 


1.2 


0502 


1202 


010 


4601 


01 


11 


6 


1.2 


0303 


09 


010 


4601 


01 


24 


3 


0.6 


0303 


09 


010 


4601 


01 


0207 


21 


4.1 


0502 


1202 


010 


4601 


01 


0207 


10 


2.0 


0502 


1401 


010 


4601 


01 


0207 


4 


0.8 


0502 


1602 


010 


4601 


01 


0207 


2 


0.4 


0303 


07 


017 


57 


0602 


01 


4 


0.8 


02 


0301 


002 


5801 


0302 


33 


26 


5.1 


06 (04-07) 


1302 


002 


5801 


0302 


33 


3 


0.6 


02 


0301 


002 


5801 


0302 


11 


4 


0.8 



15a«15 (01. 04-07, 12, 14. 19. 20, 24. 25, 26N. 27. 32^5) 
15b'<15 (02. 08, 11, 15, 28. 30) 



served MICA- HLA-B - HLA-Cw haplotypes. The three most com- 
mon three-locus haplotypes were MICA*010 - B*4601 - Cw*01 
(HF=14.3%), MICA*02 -B*5801 - Cw*0302(HF=8.2%) and 
MICA*019 - 8*15(02,08,11,15,28,30) - Cw*08 (HF=7.1%). The two 
most common six-locus haplotypes were DQB1*02 - DRB1*0301 - 
MICA*002 - B*5801 - Cw*0302 - A*33 (5.1%) and DQB1*0303 - 
DRBP09 -MICA*010 - B*4601 - Cw*01 - A*0207 (4.1%). Further, 
other than those haplotypes, as many as 38 distinct thi*ee-locus and 
28 six-locus haplotypes were recognized in more than two NET 
individuals, indicating generation of a large variety of allelic combi- 
nations among the HLA loci in NET. 



Discussion 

The MICA gene shows an extensive degree of genetic polymoi-phism. 
which is slightly more concentrated in the ot2 domain, resembling 
classical HLA class I molecules. However, in contrast to classical class 
I molecules MICA polymorphism is equally diverse within the ctS do- 
main. Further, sequence variation at each polymorphic position in the 
MICA gene has been limited to bialleiic substitution in contrast to 
many multi-allelic ones in the classical MHC class I genes. The struc- 
ture of MICA molecules revealed by X-ray crystallography confirms 
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the general configuration of the molecules with an MHC class I fold 

(20, 21). However, the MICA polymorphic residues are positioned on 
the outer edge of an antigen-binding cleft, apparently bordering an 
invariant ligand-binding site, unlike MHC class I molecules. There- 
fore, actual role of polymorphism in MICA molecules has yet to be 
determined. 

In order to investigate genetic polymorphism in the extracellular 
domains of the MICA gene in the NET population, SBT was adopted 
in this study. Consequently, MICA*008 was the most common allele, 
but its allele frequency showed a significant decrease as compared 
with those of the Japanese and Caucasian populations (see Table 2). 
The MICA*008 allele is highly dominating in Caucasians (AF=55%), 
being associated with B*0801, B*0702, 3*4402 and B*4001 (19). In 
Caucasian, its predominant abundance and representation on 
multiple haplotypes could indicate that MICA*008 is the ancesti'al al- 
lele with subsequent recombinations and point mutations producing 
further allelic divei-sity (19). In the NET population, MICA*008 (AF= 
21.4%) is mostly associated with unique HLA-B alleles (B*27 and 
B*39/6701) different from that in the Caucasian population. 
MICA*008 in the NET population carries the A5.1 allele in the TM 
(transmembrane) domain (unpublished data), as similarly observed in 
other populations (16, 22), which possibly produces a truncated MICA 
protein because of a frame-shift mutation due to a one-base insertion 
in exon 5. 

MICA*010 and *019, which are different from M1CA*008 by 
double- and single-nucleotide substitutions, are also prevalent with 
AFs of 18.2% and 15.3% in NET. respectively. These two MICA 
alleles might have been generated from MICA*008 and then con- 
served in the Thai population. However, they were associated with 
TM (A5) and HLA-B alleles (B*4601 and B*15) different from 
MICA*008 in this population. Additionally, MICA*010 carries a nu- 
cleotide substitution at position 17 in exon 2, which results in a 
proline instead of an arginine by substitution at amino acid position 
6 in the al domain. This substitution has been recently established 
to abolish cell surface expression of the MICA*010 allele by 
blocking a P-sheet hydrogen bond with the histidine carbonyl group 
at position 27 on the P2 strand. Moreover, this substitution is incom- 
patible with P-sheet secondary structure, thus likely interfering 
with protein folding (20, 23). Thus far, we have recognized six indi- 
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Genetic susceptibility to carbamazepine-induced cutaneous 
adverse drug reactions 

Shuen-lu Hung^'\ Wen-Hung Chung^'^'\ Shiou-Hwa Jee"^, Wen-Chieh Chen^ 
Yun-Ting Chang^, Woan-Ruoh Lee^, Shu-Ling Hu^, Meng-Tse Wu', 
Gwo-Shing Chen^, Tak-Wah Wong^, Pa-Fan Hsiao', Wei-Hsuan Chen^ 
Han-Yu Shih^ Wu-Hsiang Fang^ Chun-Yu Wei^ Yi-Hui Lou^ Yau-Li Huang*", 
Juei-Jueng Lin"^ and Yuan-Tsong Chen^'" 



The anticonvulsant carbamazepine (CBZ) frequently 
causes cutaneous adverse drug reactions (cADRs), 
including maculopapular eruption (MPE), hypersensitivity 
syndrome (HSS), Stevens-Johnson syndrome (SJS) and 
toxic epidermal necrolysis (TEN). We reported that 
SJS/TEN caused by CBZ is strongly associated with the 
HLA-B*1502 gene in Han Chinese. Here, we extended our 
genetic study to different types of CBZ-cADRs (91 patients, 
Including 60 patients with SJS/TEN, 13 patients with 
hypersensitivity syndrome and 18 with maculopapular 
exanthema versus 144 tolerant controls). We used 
MALDI-TOF mass spectrometry to screen the genetic 
association of 278 single nucleotide polymorphisms 
(SNPs), which cover the major histocompatibility complex 
(MHO region, tumor necrosis factor-alpha, heat shock 
protein and CBZ-metabolic enzymes, including CYP3A4, 
2B6, 2C8, 2C9. 1A2 and epoxide hydrolase 1. In addition, 
we genotyped 20 microsatellites in the MHC region and 
performed HLA-typing to construct the recombinant map. 
We narrowed the susceptibility locus for CBZ-SJS/TEN to 
within 86 kb flanking the HLA-B gene on the extended 
B*1502 haplotype, and confirmed the association of 
B*1502 with SJS/TEN [Pc = 1.6 x 10"^\ odds ratio 
(OR) =1357; 95% confidence interval (CI) = 1 93.4-8838.3]. 
By contrast to CBZ-SJS/TEN, HLA-B* 1502 association was 
not observed in the MPE or HSS groups: MPE was 
associated with SNPs in the HLA-E region and a nearby 
allele, HLA'A*3101 (Pc=2.2 x 10"^ OR = 17.5; 95% 
Cl=4.6-66.5), and HSS with SNPs in the motilin gene 
(Pc=0.0064, OR = 7.11; 95% 01 = 3.1-16.5) located terminal 
to the MHC class II genes. No SNPs in genes involved in 
CBZ metabolism were associated with CBZ-induced 
cADRs. Our data suggest that HLA'B*1502 could 
contribute to the pathogenesis of CBZ-SJS/TEN, and that 



Introduction 

Although carbamazepine (CBZ) is a commonly prescribed 
first-line anticonvulsant, it is also a relatively common 
cause of cutaneous adverse drug reactions (cADRs) [1]. 
The frequency of CBZ-induced cADRs is between 1 in 
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1000 and 1 in 10000 new exposures to the drug in 
Caucasians [2,3]. The cADRs range from mild maculo- 
papular eruption (MPE), with increasing severity, to 
hypersensitivity syndrome (HSS), Stevens-Johnson syn- 
drome (SJS) and toxic epidermal necrolysis (TEN). The 
mortality rate of the latter can be as high as 40% [4]. 
MPE is characterized by cutaneous fine pink macules and 
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papules, lesions which usually fade within 1-2 weeks 
following cessation of drug treatment. HSS is character- 
ized by multi-organ involvement (e.g. hepatitis and 
nephritis) accompanied by systemic manifestations (e.g. 
fever, arthragias, cosinophilia and lymphadenopathy) in 
addition to skin rashes [4]. Skin manifestations of HSS 
may vary from MPE to exfoliative dermatitis [5]. HSS is 
also called drug rash with cosinophilia and systemic 
symptoms. SJS and TEN are characterized by a rapidly 
developing blistering exanthema of purpuric macules and 
target-like lesions accompanied by mucosal involvement 
and skin detachment [6]. According to the clinical 
morphology, SJS/TEN belong to the group of bullous 
cADRs, whereas MPE and HSS are non-bullous reactions. 

The mechanism by which CBZ causes cADRs is not well 
understood. Potential defects in the enzymes responsible 
for bioactivation and detoxification of CBZ have been 
proposed [1]. Carbamazepine is bioactivated by hepatic 
cytochrome P450 enzymes, mainly CYP3A4, CYP2B6 and 
CYP2C8, which generate various potentially reactive 
metabolites, such as CBZ 10,1 1-epoxide, 3-hydroxy 
CBZ, 2-hydroxy CBZ, and CBZ 2,3-epoxide [7,8]. Most 
of the reactive epoxides are detoxified to non-toxic 
dihydrodiols by liver microsomal epoxide hydrolase 1 
(EPHXl) or to glutathione conjugates by glutathione 
transferase [9,10]. Previous attempts to identify a genetic 
defect altering the structure or function of epoxide 
hydrolase 1 in individuals susceptible to CBZ-induced 
cADRs, however, were not successful [11,12]. 

It has also been proposed that immune reactions are 
involved in the cADRs triggered by CBZ because 
infiltrating inflammatory cells can be detected in the 
skin lesions, and rechallenging with the same drug 
typically shortens the incubation period, resulting in 
more severe manifestations [13,14]. CD4 + T cells are 
the major cell type found in the skin lesions of MPE 
and HSS [15], whereas CDS + T-cell-mediated 
cytotoxic responses appear to be the major event in 
SJS/TEN [16]. There is also evidence supporting the 
view that cADRs involve major histocompatibility com- 
plex (MHC)-dependent presentation of its metabolites 
forTcell activation [16,17]. Naisbitt ^/ [15] reported 
that CBZ might be presented by MHC class II expressed 
on the surface of antigen-presenting cells to the T-cell 
receptor of CD4 + T cells of patients with CBZ-induced 
HSS [15]. 

We recently reported that, in Han Chinese, the HLA- 
8*1502 g^nG is strongly associated with CBZ-induced SJS 
and TEN [18]. In the present study, we sought to extend 
our genetic association study to other types of cADRs 
induced by CBZ using a comprehensive genetic screen of 
markers in the MHC region, as well as genes encoding 
CBZ metabolic enzymes. In addition, we performed fme 



mapping to further investigate the susceptibility locus for 
CBZ-induced SJS/TEN. 

Materials and methods 
study population 

From 1997 to 2004, 91 individuals who fulfilled diagnostic 
criteria for CBZ-induced MPE, HSS, SJS or TEN agreed 
to participate in the study. Eighty-eight of the 91 patients 
were recruited from Chang Gung Memorial Hospital 
Health System, Veterans General Hospital, National 
Taiwan University Hospital, Cathay General Hospital 
and Chushang Show-Chwan Hospital, Taiwan. We 
recruited patients primarily from inpatients who had 
more severe diseases. The remaining three patients were 
referred by physicians from Hong Kong (two patients) 
and the USA (one patient). The latter three all had CBZ- 
SJS. All patients (including the three from Hong Kong 
and the USA) and controls (see below) were Han 
Chinese or Chinese descendants. 

The 91 patients included 44 patients with CBZ-SJS/TEN 
who we previously reported [18], and 16 additional 
patients with SJS/TEN, as well as 13 patients with 
hypersensitivity syndrome and 18 with maculopapular 
exanthema. All patients were assessed by two dermatol- 
ogists who reviewed photographs, pathological slides and 
medical records. Diagnostic criteria for SJS/TEN were 
based on the clinical morphology defined by Roujeau [6]. 
We defined *SJS' as a skin detachment of less than 10% of 
total body-surface area, 'overlap of SJS and TEN* as skin 
detachment of 10-30%, and 'TEN' as skin detachment 
greater than 30%. The criteria for HSS in this study were 
skin rash, plus two of the following symptoms: fever, 
lymphadenopathy and haematologic abnormalities (e.g. 
cosinophilia, atypical lymphocytosis) with involvement of 
at least one internal organ (e.g. hepatitis, pneumonitis, 
myocarditis, pericarditis, nephritis) [2,5,19]. SJS and 
TEN are bullous cADRs and are considered to be variants 
of the same disease. Therefore, they were analysed 
together as a group. MPE and HSS are non-bullous and it 
is less clear whether they are from a single disease 
spectrum. They were therefore analysed as a group and 
also analysed separately as two entities. In all enrolled 
cases, CBZ was regarded as the offending drug if the 
onset of cADRs symptoms occurred within the first 2 
months of exposure and the symptoms resolved upon 
withdrawal of the drug. Patients with an absence of 
symptoms after re-exposure to CBZ were excluded. 

The control group was the 144 consecutive patients who 
received CBZ for at least 3 months without evidence of 
adverse reactions. These tolerant patients were recruited 
from the Neurology Clinic of the same regional hospitals 
where cADRs patients were recruited. In addition, 93 
healthy subjects were randomly selected from a biobank 
under a nationwide population study, in which 3312 Han 
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Chinese descendants were recruited based on the 
geographical distribution across Taiwan. There was no 
self-report of adverse events in any of these 93 subjects. 
The study was approved by the institutional review 
board, and informed consent was obtained from ait of the 
participants. 

DNA isolation and genotyping 

Genomic DNA was isolated using the PUREGENE DNA 
purification system (Gentra Systems, Minneapolis, Min- 
nesota, USA). 

Short tandem repeat polymorphism (5TRP) genotyping 

Twenty highly polymorphic microsatellite markers 
located in the MHC region were selected from the 
NCBI database (i.e. D6S258, D6S2972, D6S510, 
D6S265, D6S388, D6S2814, HLAC-GA, HLABC-CA, 
MIB, MICA, TNFd, BAT2-CA, D6S273, D6S1615, 
DQCAR, G51152. D6S2414, D6S1867, D6S1560 and 
D6S1583). The average heterozygosity of markers was 
0.702 with an estimated 230 kb of spacing. Primers were 
designed based on oligonucleotide sequences reported 
within the database. Polymerase chain reaction (PGR) for 
genotyping was performed in a 5-fil volume containing 
10 ng of genomic DNA and 0.33 |iM of each primer by 
using GeneAmp 9700 thermocyclers (Applied Biosys- 
tems, Foster City, California, USA). Up to six products of 
appropriate size and carrying a fluorescent label were 
pooled for capillary gel electrophoresis. The size of 
polymorphic amplicons was determined by ABI 3730 
DNA sequencer (Applied Biosystems) using the LIZ500 
size standard as an internal size standard. Allele sizing 
was calculated using the GENMAPPER program version 
3.0 (Applied Biosystems). Allele calling and binning were 
performed using the SAS program (SAS Institute, Gary, 
North Carolina, USA). Three CEPH control individuals 
(1331-01, 1331-02, 1347-2) and H2O were included in all 
genotyping experiments for quality control purposes. 

Single nucleotide polymorphisms (SNP) genotyping on 
IMHC region and metabolic enzymes for CBZ 

A total of 379 SNPs were selected from the SNP database 
(dbSNP: build 123) of the National Center for Biotech- 
nology Information for genotyping [20]. These included 
220 SNPs from 4Mb of the MHC region on chromosome 
6p21.3 and 159 SNPs selected from genes encoding drug 
metabolizing enzymes. The 220 SNPs selected in the 
MHC region included 201 SNPs reported by Walsh e^a/. 
[21], and additionally included SNPs of tumor necrosis 
factor (TNF) -alpha and the heat shock proteins. We 
especially included rsl800629, TNF-308A allele, which is 
reported to be associated with CBZ-HSS [22], and 
rs2227956 of HspAlL (heat shock 70 kDa protein l-like), 
which is reported to be associated with abacavir 
hypersensitivity [23]. The average space of selected 
SNPs in the MHC region was approximately 20 kb. 
Another 159 SNP set contained genes potentially 



involved in bioactivation or detoxification of CBZ 
metabolism, including CYP3A4, CYP2B6, CYP2C8, 
CYP2C9, CYP1A2 and EPHXl. 

SNP genotyping was performed using high-throughput 
MALDI-TOF mass spectrometry. Briefly, primers and 
probes were designed using SpectroDESIGNER software 
(Sequenom, San Diego, California, USA). Multiplex 
PCRs were performed, and unincorporated dNTPs were 
dephosphorylated using shrimp alkaline phosphatase 
(Hoffman-LaRoche, Basel, Switzerland) followed by 
primer extension. The purified primer extension reaction 
was spotted onto a 384-element silicon chip (Spectro- 
CHIP, Sequenom), analysed in the Bruker Biflex III 
MALDI-TOF SpectroREADER mass spectrometer (Se- 
quenom) and the resulting spectra processed with 
SpectroTYPER (Sequenom). 

Of the 379 SNPs selected from the database, 101 SNPs 
were excluded from further analysis due to the non- 
polymorphism, low successful rate (less than 0.7), and/or 
departure from Hardy-Weinberg equilibrium iP< 0.001) 
in our control group. The remaining 278 SNPs, including 
88 SNPs in metabolic pathway for CBZ and 190 SNPs in 
the MHC region, were used for the statistical analyses 
and the P-values of association were adjusted by using 
Bonferroni correction. 

HLA genotyping 

HLA alleles A, By C and DRBl were determined by 
sequence-specific oligonucleotide reverse lineblots (DY- 
NAL Biotech Ltd, Bromborough, UK) (24). Potential 
ambiguities were resolved by sequencing-based typing 
[25]. The primers used for PGR amplification of DNA 
fragments of HLA-A, B and C genes were: (i) Binl-TA- 
M13F, Binl-CG-M13Fand Bin3-M13R for T/Z^-^ alleles; 
(ii) Ainl-A-M13F, Ainl-G-M13F, Ainl-T-M13Fand Ain3- 
62-M13R for HLA-A alleles; and (iii) 5Clnl-61 and 
3BCln3-12 for //Z^- 6* alleles [26]. The DRBl alleles were 
separated into eight groups: DRl, DR2, DR3/11/6, DR4, 
DR7, DR8/12, DR9, DRIO for sequencing-based typing 
as described [27]. Genotyping for MHC class I chain- 
related gene A (MICA) was carried out by DNA 
sequencing on its exons 2, 3, 4 and 5 [25]. The 
primers used for PGR amplification of the DNA fragment 
of the MICA gene were: forward primer: 5'- 
CGTTCTTGTCCCTTTGCCCGTGTGC-3' and re- 
verse primer: 5'-GATGCTGCCCCATrCCCTTCCCAA-3'. 
The sequence data were analysed by SeqScape v2 
(Applied Biosystems). 

Statistical analysis 

Comparisons of allele or genotype frequencies between 
groups were performed using Fisher's exact tests. All 
Avaiues were two-tailed. P < 0.05 was considered to be 
statistically significant. An allelic association screen was 
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carried out by the Cochran-Armitage Trend test for each 
SNP or STRP marker [28]. To test the association of 
haplotype frequencies, the composite haplotype method 
and haplotype trend regression of Helix Tree software 
version 3.0.0 were used. The linkage disequilibrium 
between two loci was analysed by using Expectation/ 
Maximization (EM) method provided by Helix Tree 
software version 3.0.0 (Golden Helix Inc., Bozeman, 
Minnesota, USA). To achieve sufficient power to identify 
loci associated with different clinical manifestations, the 
corrected P {P^) values were adjusted by using Bonferro- 
ni's correction for multiple comparisons (278 for SNP 
assays, 17 for HLA-A, 40 for HLA-B, 19 for HLA-C, 30 for 
HLA-DRBl, 25 for MICA, and 20 for HLABC-CA). Odds 
ratios (ORs) were calculated using Haldane's modifica- 
tion, which adds 0.5 to all cells to accommodate possible 
zero counts [29]. 

Results 

Characteristics of patients and controls 

Based on the clinicopathologic features of cADRs, our 
patients could be divided into bullous and non-buUous 
groups: (i) bullous: CBZ-SJS/TEN {n = 60; 54 with SJS, 
five with overiapping SJS/TEN, one with TEN) and (ii) 
non-bullous: GBZ-MPE/HSS (/7 = 31; 18 with MPE, 13 
with HSS). The non-bullous group was further separated 
into two subgroups (MPE and HSS) for data analysis. 
Clinical manifestations and demographic variables of the 
91 patients and 144 tolerant controls are summarized in 
Table 1. The onset of symptoms for all patients with 
cADRs occurred within the first 2 months of CBZ 
exposure. The mean duration of CBZ exposure was longer 
in HSS patients than in SJS/TEN. Three patients with 
SJS/TEN had a second attack within 2 days of re- 
exposure; one of them developed TEN and died during 
the second attack. All 60 patients with SJS/TEN had 
widespread purpuric rash with blisters and mucosal 
involvement. In addition to maculopapular exanthema, 



all 13 patients with HSS had symptoms of multi-organ or 
systemic involvement. 

One hundred and forty-four patients who had been on 
CBZ for at least 3 months (mean 75.03 months) with no 
self-report of adverse events were used as tolerant 
controls. Their demographic variables are shown in Table 
1. Note that the CBZ-tolerant controls received higher 
doses of CBZ (mean dosage 793.8 mg/day) and yet no 
adverse drug reactions were observed. 

Fine mapping of the major histocompatibility 
complex region for genetic susceptibility to 
CBZ-SJS/TEN 

To confirm our previous observation of the genetic 
association between HLA-B*1502 and CBZ-SJS/TEN, 
and to further define the suscepubility region, we 
genotyped 20 STRP markers and 190 SNPs in the 
MHC region. Nine STRP markers located between 
D6S2814 and D6S2414 (physical map 30.8-32.9 Mb, 
chromosome 6) showed significant association 
(/'-values < 10"^'; Fig. 1). In particular, HLABC-CA near 
HLA'B showed the strongest association {P — 3.4 x 10" 
Fig. 1), followed by D6S1615, a marker closed to HLA- 
DRBl locus. 

Thirty-seven SNPs showed /^-values of less than 0.001 
(Fig. 2a). Among them, seven SNPs located between 
HLA-DRA and HIA-C (physical position 31.3-32.5 Mb, 
chromosome 6) showed strong associations (P < 10"^^) 
(Fig. 2a). The three most important SNPs were 
rs3130690, rs2848716 and rs750332, located near HLA- 
B, MICA and BATS genes (Table 2). The rs3 130690, an 
intergenic SNP with 36 kb telomeric to the HLA-B locus, 
demonstrated the strongest association with CBZ-in- 
duced SJS/TEN (/>c= 1.29x 10"^^^) (Table 2). The T 
allele of the rs3 130690 SNP was present in 95% (57/60) 
of CBZ-SJS/TEN patients, but only in 6.9% (10/144) of 



Table 1 Demographic variables, dosage and duration of carbamazepine (CBZ) exposure, and clinical characteristics in CBZ-induced 
cutaneous adverse reactions and CBZ-tolerant controls 

Bullous cADRs Non-bullous cADRs (n = 3 1 ) 

SJS/TEN (n = 60) HSS (n « 1 3) MPE (n « 1 8) Tolerant controls (n = 1 44) 



Sex, n (%) 



Male 


33 (55) 


9 


Female 


27 (45) 


4 


Age (years), mean (range) 


43.4 (5-80) 


51.5 


CBZ exposure, mean (range) 






Dosage (mg/day) 


332.6 (100-600) 


420 


Duration 


15.1 (2-49) days 


32.7 


Cutaneous features, n (%) 






Blister or epidermal detachment 


60 (100) 


0 


Mucosal erosions 


60 (100) 


3 


General and laboratory, n (%) 






High fever OSS-B^C) 


41 (68.3) 


12 


Eosinophil count > 1 000/mm^ 


4 (6.7) 


10 


Atypical lymphocytosis 


8 (13.3) 


7 


Abnormal liver function 


7 (11.7) 


11 


Abnormal renal function 


1 (1.7) 


2 



(69.2) 


8 (44.4) 


75 (52.1) 


(30.8) 


10 (55.6) 


69 (479) 


(6-83) 


45,9 (7-84) 


35.7 (5-79) 


(300-600) 


423.1 (200-800) 


793.8 (100-1500) 


(15-54) days 


22.4 (7-55) days 


75.03 (3-2875) months 


(0) 


0(0) 


Not observed 


(23) 


0(0) 


Not observed 


(92.3) 


11 (61.1) 


Not observed 


(76.9) 


4 (22.2) 


Not obsen/ed 


(53.8) 


3 (1 6.7) 


Not observed 


(84.6) 


0(0) 


Not observed 


(15.4) 


0(0) 


Not observed 



SJS/TEN, Stevens-Johnson syndrome/toxic epidennal necrolysis; HSS, hypeisensitivity ^ndrome; MPE, maculopapular exanthema. 
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Association screen of 20 short tanden^ repeat polymorphism (STRP) markers in the major histocompatibility complex (MHC) region with 
carbamazepine (CBZ)-induced cutaneous adverse daig reactions. On the ;r-axis, 20 STRP markers in the MHC region are ordered by their physical 
positions (29.9-33.9 Mb) on chromosome 6p21.3. On the /-axis, the -log io /^values were calculated by comparison of the allele frequencies 
between the patients and tolerant group using the Cochran-Anrtitage exact trend test. Genotyping data of four groups of patients are presented: 
diamond symbols: CBZ-induced Stevens-Johnson syndrome (SJS)/toxic epidermal necrolysis (TEN); square symbols: CBZ-induced maculopapular 
eruption (MPE)/hypersensitivity syndrome (HSS); triangle symbols: CBZ-induced HSS; and circle symbols: CBZ-induced MPE. 



tolerant controls [OR = 254.6; 95% confidence interval 

(CI) =70.4^-901.9] for TT/TG genotype versus GG 
genotype) (Fig. 3). To verify the results of our association 
screen, we also genotyped 93 healthy subjects randomly 
selected from the general population [30]. A strong 
association between SNPs and CBZ-induced SJS/TEN 
was replicated in the MHC region, especially rs31 30690, 
nearby the HLA-B locus (data not shown). 

We further genotyped the individual HLA-A, 5, C, DRBl 
and MICA alleles. The HLA-B*I502 allele was present in 
98.3% (59/60) of CBZ-SJS/TEN patients, whereas only 
4.2% (6/144) of the tolerant controls were positive for the 
allele {Pc=l.6x lO'^^^ OR = 1357; 95% CI =193.4- 
8838.3) (Table 3 and Fig. 3). The only CBZ-SJS patient 
who did not have the B*1502 allele, had instead, another 
HLA-B15 allele: HLA-B*I558, In addition to the positive 
association, we also found that NLA'B*400J was negatively 
associated with CBZ-SJS/TEN {P^^2.6x 10"^ 0R = 
0.16; 95% CI = 0.1^.4) (Table 3). The NLA-Cw*080J and 
MICA*0I9 alleles flanking the HLA-B'^1502 showed strong 
linkage disequilibrium with HLA'B*I502, and were present 
in 93.3% (56/60) and 95% (57/60) of SJS/TEN patients, 
respectively. An extended B*1502 haplotype formed by 
polymorphic alleles (A*1101-Cw*0801-HLABC-CA*119- 
rs3130690T-B*1502-MIGA*019-DRBl*1202) had a strong 
association with CBZ-induced SJS/TEN. The recombinant 
map of Cw*0801-HLABC-GA*119-rs3130690T-B*1502- 
MICA*019 defined the susceptible region within 86 kb 
(i.e. between the T allele of rs3 130690 and MICA* 01 9) 



flanking the B*I502 gene in the 4Mb MHC region 
(Fig. 3). Within this 86-kb region, HLA-B is the only 
known gene. Taken together, the data suggested that one 
or more alleles in the vicinity of the HLA-B locus, 
particularly B*I502 itself, participate in the pathogenesis of 
CBZ-induced SJS/TEN. 

Association screen for candidate gene SNPs with 
CBZ-iVIPE/HSS 

When MPE and HSS were grouped together as non- 
bullous cADRs, no STRP markers showed association, 
and only six SNPs had P-values less than 0.01 (Figs 1 and 
2b). The three most significant SNPs are listed in 
Table 2, including rsl264511 located near HLA-E, 
rsl042389 located on the 3'-untranslation region of 
CYP2B6, and rs2894342 located in the promoter region 
of the motilin (MLN) gene (Table 2). However, the 
P-values became non-significant after correcting for 
multiple testing (278 SNP assays) was performed. 

In HLA-typing, HLA-A*3I0I showed an association with 
MPE/HSS (/^c = 0.0021) (Table 3). HU-A*3I0I was 
present in 25.8% (8/31) of patients with MPE/HSS, but 
only in 2.8% (4/144) of tolerant controls (OR =12.17; 
95% CI = 3.6-41.2). 

Association between MLN polymorphisms and HSS 

When genetic analysis was performed separately for HSS 
and MPE, one STRP marker, D6S2414 (physical map: 
32.9 Mb, chromosome 6) showed an association with 



Copyright © Lippincott Williams & Wilkins. Unauthorized reproduction of this article is prohibited. 



302 Pharmacogenetics and Genomics 2006, Vol 1 6 No 4 



CA 

c 
o 

(0 

9 



0) 

in 



1 i 

to ^ 



73 

c 2 

O) t: _ 
tf) M .9 

^ o .tr 
o g a 

! 5 



CO ^ 
*- CD 



I I I I I I I I I I 
CO 1-^ ^' CO r-" *^ r-' ^* 



o o o 

X X XZZZoZZZZZ 

o o to d 

CM -^t 

*^ oi ID 



CM CM 

I I "incD cocO'-cocot^ 

OOOC0<»CNO<t<-00 00P^<» 

^_;ooooooooo 

tOflO^OOOOOOOOO 

!2!2Sodd6ddddd 
^ « 



m(oa>ooa>ooa>cMoa)cooo 
■^^coooaoocscDO^ro 

dddddddddddd 



CO»-C0tJ-«JQ0CMCOCOCOCM00 

oot»oc^t-m«)cMcNcr>t<-e*. 
(otomcMincM^^^co^oi 
d d dddododddd 



v-0Oe>O<<C!J0OHO 



c ojq; 



CD m LU 

3 < <o ^ S 2 2 5 3 3 4- - 
I y ^ ^ £ 1 1 £ i ^ 3 ■ 



CO 



»- 

O) Oi 

CO lO 

a> o> 

CO ^ 

CO CO 

«9 tt> 



OS ■»- 
CM in 

o O) 

in ^ 

CM 
lO 

^ o 

CO CO 



S CM 



CM lO ^ T- »- CO 

CM in in 00 

CO O) (O <7) (D 

CM in ^ ^ in -St 

CO CO CO CM U} O) 

00 CO 00 in in CD 

CO T- CO o o f- 

C3 CO W CO CO CO 

<0 CO O <D (D CO 



OcOrvi^CBCMCMOCM^Oin 

cDf-Sincococooo^inin^ 

ooogTrcM-ia-^tininTj-cnS 
co'*Jj;<£>'«3-c»c3)t>.<j>cDin5j 

T-CDH2c>IOOOflOOCOCMOS 
COC>IVX-t--»-CMCMCMCMT-^<S 



w 
I 

— . 
UJ 



CO 
CO 

I 



CL 



E 

8 

i 

E 



(0 • 



fl> e 
•2, §; 

(0 

p CO 

IS 

li 

or § 



S E 

"i- Q. 

^2 



1 j5 
o 

V o 
o s 
c o 

0) < 
C E 

2 o 



§co 

o»ln 
g> o 
0) d 
5 A 

er 

o § 

9M 



: CO 
z 



CBZ-induced HSS {P=6.2S x 10-^ Fig. 1). In addition, 
several SNPs located in the MHC class II terminal region, 
around the motilin.gene (physical map: 33.8 Mb), showed 
associations with HSS (Table 2 and Fig. 2c). The most 
significant SNP rs2894342 was located in the promoter of 
the MLN gene, and its A allele had an increased risk for 
CBZ-HSS (Pc = 0.0064, OR =7.11; 95% CI = 3.1-16.5) 
(Table 2), The next most significant SNP was a 
nonsynonymous SNP (rs2O75800 of HSPAIL, heat shock 
70 kDa protein 1-like) involving an A to G transition 
which leads to an amino acid change of lysine to 
glutamate at residue 602. The third most significant 
SNP was rs2395402, which was located in the intron of 
the LEM domain containing 2 (LEMD2) gene, 21 kb 
telomeric to rs2894342 of the MLN gene (Table 2). 
However, the P-values for the latter two SNPs became 
non-significant after Bonferroni's correction. 

None of the HLA-A, B, 6* and DRBl alleles had significant 
associations with CBZ-induced HSS (Table 3). 

Association between HLA'A*3101 and MPE 

STRP or SNP markers in the HLA-B or MLN region on 
chromosome 6p21.3 showed no association with MPE 
(Figs 1 and 2d). Instead, two SNPs near the HLA-E 
region (physical position 30.5 Mb) and one SNP in the 
5'-untranslation region of the leukocyte specific transcript 1 
(LSTl) gene showed associations with MPE 
(Fig. 2d and Table 2, P = 0.00078-0.0075). The two 
SNPs in the HLA-E region were: rsl26451 1, an intergenic 
SNP located near HLA-E and rs 10595 10 located in HLA- 
E gene (Table 2). HLA-genotyping further revealed that 
the HLA'A*3J0I allele, approximately 530 kb telomeric to 
HLA-E, was present in 33.3% (6/18) of MPE patients but 
only in 2.8% (4/144) in the tolerant control group 
{Pc=Z2x 10'\ OR =17.5; 95% 01 = 4.6-66.5) 
(Table 3). 

Association screen for SNPs In drug metabolizing 
enzymes 

Genotyping of 88 informative SNPs in genes that are 
potentially involved in metabolic pathways for CBZ, 
including CYP3A4, CYP1A2, CYP2C8, CYP2B6, CYP2C9 
and EPHXl, in various types of CBZ-induced cADRs, 
showed that only one SNP, rsl042389, in the 
3'-untranslation region of CYP2B6, had a weak association 
with CBZ-MPE/HSS (P = 0.0016, Table 2 and Fig. 2b). 
The P-value became non-significant after Bonferroni's 
correction, suggesting that the genetic polymorphisms of 
genes encoding drug metabolizing enzymes have no 
significant association with CBZ-induced cADRs, regard- 
less of the type of cADRs, 

Discussion 

This study confirmed our previous observation that HLA- 
8*1502 is strongly associated with CBZ-SJS/TEN in Han 



Copyright © Lippincott Williams & Wilkins. Unauthorized reproduction of this article is prohibited. 



Carbamazepine-induced adverse drug reactions Hung et aL 303 



Fig. 2 




Screening of 278 candidate single nucleotide poiynriorphisms (SNPs) for association with (a) carbamazepine (CBZ)-induced Stevens-Johnson 
syndrome/toxic epidermal necrolysis, (b) CBZ-induced maculopapular eruption (MPE)/hypersensitivity syndrome (HSS), (c) CBZ-induced HSS and 
(d) CBZ-induced MPE On the r-axis, 278 SNPs are ordered by their chromosomal positions, including 58 SNPs of epoxide hydrolase 1 (EPHX1 ) on 
chromosome 1 , 1 90 SNPs in the major histocompatibility complex (MHC) region on chromosome 6, and 30 SNPs of CYP3A4, 2C8, 2C9, 1 A2 and 
2B6. On the /-axis, the -log^o P-values were calculated by comparison of the allele (diamond symbols) or genotype (triangle symbols) frequencies 
between the patients and tolerant group using the Cochran-Armitage exact trend test. 



Chinese [18]. We increased the number of study patients 
from the original 44 to 60 and included patients from a 
wide geographical distribution (Taiwan, Hong Kong, 
China and the USA). However, all patients enrolled were 
Chinese or Chinese descendants. The fact that the HLA- 
B* 1502 rWgIc was present in low frequency in Caucasians 
may explain the low incidence of CBZ-SJS/TEN in 
Caucasians [31,32]. 

It is interesting to note that the only CBZ-SJS/TEN 
patient who did not have the HLA'B*J502 gene had 



HLA'B*I558y another HLA-B15 variant. This variant 

allele was present in low frequency in our population 
(approximately 0.9%) and was not detected in any of the 
144 tolerant patients. The detection of B*J558 in one of 
our patients may imply that B*I558 shares a similar 
structural feature with B*I502 for triggering the immune 
reaction of SJS caused by CBZ. 

Fine recombinant genetic mapping on the extended HLA- 
B''1502 haplotype (A*l 101-Cw*0801-HLABC-GA*l 19-rs 
3130690T-B*1502-MICA*019-DRB1*1202) further nan-owed 
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Fig. 3 
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Susceptible region 

Susceptible region for carbamazepine (CBZ)-induced Stevens-Johnson syndrome (SJS)/toxic epidermal necrolysis (TEN) on the extended HLA- 
8*1502 haplotype. Schematic map of 4 Mb major histocompatibility complex region on chromosome 6p21 .3 shows the relative positions of HLA-A, 
B, C, DRB1 and MICA genes, STRP marker (HLABC-CA), and SNP marker (rs31 30690). The percentage of presence of the specific allele in CBZ- 
SJS/TEN patients and tolerant-controls is shown. The susceptible region should be within 86 kb (I.e. between T allele of rs3 130690 and 
MICA*019). Within this region, HLA-B* 1502 is the only known gene and shows the strongest association with CBZ-induced SJS/TEN. NS, Not 
significant (P>0.05). 



Table 3 Association of HLA alleles with carbamazepine (CBZ)-induced cutaneous adverse reactions 



Tolerant controls SJS/TEN (n= 60) MPE/HSS (n=31) HSS (n=13) MPE(n=18) 

(n=144) 



HLA-A 



*1101 


67 


42 [NS; 2.68 (1.4-5.1)] 


12 [NS; 0.73 (0.3-1.6)] 


6 [NS; 0.99 (0.3-3)] 


6 [NS; 0.58 (0.2-1.6)1 


*2402 


41 


5 [0.026; 0.23 (0.1-0.6)1 


9 [NS; 1.03 (0,5-2.4)1 


3 [NS; 0.75 (0.2-2.7)] 


6 [NS; 1.26 (0.5-3.5)1 


*3101 


4 


1 [NS;0.59 (0.1-4.1)1 


8 [0.0021; 12.17 (3.6-41.2)1 


2 [NS; 6.36 (1.2-33.9)1 


6 [2.2 X 10"^; 17.5 (4.6-66.5)1 


HLA-B 












*1502 


6 


59 [1.6 X 10"*^ 1357 
(193.4-8838.3)] 


1 [NS;0.77 (0,1-5.1)1 


0 [NS; 0.79 (0.1-8.8)1 


1 (NS; 1.35 (0.2-9.3)1 


*4001 


59 


6(2,6x10"*; 0.16(0.1-0.4)1 


10 [NS; 0.69 (0,3-1.5)] 


3fNS; 0.43 (0.1-1.5)) 


7 (NS; 0.92 (0.4-2.4)1 


HLA-Cw 












*0102 


54 


8(8.6x10-^ 0.26 (0.1-0.6)1 


14 [NS; 1.4 (0.6-3.0)) 


5 [NS; 1.04 (0.3-3.2)1 . 


9 [NS; 1.67 (0.6-4.3)] 


*0801 


20 


56 [7.8 X 10""; 86.8 (29.3-254.6)1 


2 INS; 0.43 (0.1-1.8)1 


0 [NS; 0.23 (0-2.3)] 


2 [NS; 0.78 (0.2-3.3)] 


HLA-DRB1 












*0405 


25 


1 [0.03; 0.08 (0.01-0.5)1 


8 [NS; 1.66(0.7-4.1)1 


1 [NS; 0.40 (0-2.5)1 


7 [NS; 3.03 (1.1-8.4)1 


*1202 


23 


41 [2.3x 10~^\ 11.4 (5.6-22.9)1 


5 [NS: 1.01 (0.4r2.8)l 


3 [NS; 1.58 (0.4-5.8)1 


2 [NS; 0.66 (0.2-2.8)] 



Data are genotype, n of positive subjects (Pc value; odds ratio (95% confidence interval)]. The association of HLA-alleles was examined by Fisher's exact test and the Pc 
values were adjusted by using Bonfenx)ni's con-ection for multiple comparisons (17 for HLA-A, 40 for HLA-B, 19 for HLA-C and 30 for HLA-DRB1). NS, Not significant 
(P>0.05). SJS/TEN, Stevens- Johnson syndrome/toxic epidermal necrolysis; MPE, maculopapular exanthema; HSS, hypersensitivity syndrome. 



down the susceptibility region for CBZ-SJS/TEN to 
within 86 kb flanking the HLA-B gene. In this region, 
HLA-B is the only gene present. Thus, it is most likely 
that B*1502 itself is the susceptibility gene for CBZ-SJS. 
Furthermore, it has been proposed that the molecular 
mechanism causing SJS/TEN involves CDS + cytotoxic 



T-cell-mediated cell death [16]. The HLA-B gene can 
elicit immune responses by presenting endogenous 
antigens to the cytotoxic T cells and it has been proposed 
that the drug and/or its metabolites may bind to the 
peptides which are then presented by the HLA-B allele 
and recognized by the specific T cells [33,34]. Our 
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cytotoxic assays also supported the view that HLA'B*I502 
is functionally involved in the cytotoxic response 
mediated by the activated T cells in CBZ-induced SJS 
(unpublished data). The recent findings of HLA'B*5101 
in abacavir hypersensitivity [35] and HLA'B*5801 in 
allopurinol-induced severe cADRs [26] further support 
the important role of the HLA-B gene in these serious 
drug-induced adverse conditions. 

SJS and TEN are cADRs characterized by bullous lesions. 
We observed a single patient who developed a bullous 
fixed drug eruption while on CBZ. Interestingly, this 
patient was also positive for HLA-B*I502 (data not 
shown). However, the present study clearly showed that 
non-bullous cADRs caused by CBZ are not linked to the 
HLA-B*I502 marker, suggesting that the genetic associa- 
tion of GBZ-cADRs is phenotype-specific. It should be 
pointed out that our previous study of allopurinol induced 
cADRs showed that both SJS/TEN and HSS were 
associated with a same HLA-B allele [26]. However, in 
this study, CBZ-HSS is not linked to the HLA-B marker 
as SJS/TEN. At the present time, we do not know why 
this discrepancy occurs. Other genetic or environmental 
factors may contribute and lead to various spectrum of 
adverse reactions caused by different drugs. 

Although there is no direct evidence that MPE and HSS 
are the same disease spectrum with differences in 
severity, HSS and MPE may be grouped together as 
non-bullous cADRs because of the similarity of the 
cutaneous manifestation [2,5,19]. However, the present 
study suggests that HSS and MPE are two disease 
entities because they are linked to different genetic 
markers: CBZ-MPE was associated with HLA'A*3101, 
and CBZ-HSS with the motilin gene polymorphisms in 
the MHC class II terminal region (Figs 1 and 2 and 
Tables 2 and 3). We did not observe the association of 
CBZ-HSS with the TNF2 allele (rs 1800629; TNF-308A 
allele; physical position: 31.6 Mb; Fig. 2c) or the ancestral 
haplotype 8.1 on the MHC, as reported previously in a 
study of CBZ-induced HSS in Caucasians [22]. This 
could be due to different study populations or different 
criteria in delineating the clinical phenotypes. In the 
present study, the SJS/TEN patients outnumber the 
patients with MPE as a reflection of the recruitment 
process because we recruited patients primarily from 
inpatients who had more severe diseases. However, in 
Han Chinese, patients with CBZ-SJS/TEN were indeed 
more frequently seen than CBZ-HSS (data not shown). 
Because of the small numbers of the patients enrolled in 
this study, further studies with a larger sample size will be 
needed to confirm this initial observation. 

Naisbitt ei al. [15] cloned T cells of patients with CBZ- 
MPE/HSS and observed that T cells may recognize CBZ 
depending on the presence of HLA class II (DR/DQ)- 



matched antigen-presenting cells. Interestingly, in our 
SNP and STRP screen and HLA typing data, we observed 
associations between CBZ-induced HSS and polymorph- 
isms markers located terminal to MHC class II region. 
However, the association of HLA'A*310J with MPE 
appears in disagreement with the observation that 
CD4 + T cells are the major cell type found in the skin 
lesions of MPE [15]. This could suggest that the gene 
associated with MPE lies in the vicinity of the HLA-A 
locus which is in linkage disequilibrium with HLA- 
A*3J0I. 

In conclusion, our data show that genetic susceptibility to 
CBZ-induced cADRs is phenotype-specific. The suscept- 
ibility gene for CBZ-induced SJS/TEN lies within an 
86 kb region flanking the HLA-B locus and HLA-B* 1502 
itself could be directly involved in the pathogenesis of 
CBZ-induced SJS/TEN. The tight association of NLA- 
B*I502 and CBZ-SJS/TEN provides a plausible basis for 
the development of such a test to identify individuals at 
risk for this potennally life-threatening condiuon caused 
by CBZ in Han Chinese, as well as for a further increased 
understanding of the pathogenesis of the clinical 
syndrome. 
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